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ABSTRACT 

Predictions of the microwave thermal emission from the interplanetary dust cloud are made using several contemporary meteoroid 
models to construct the distributions of cross-section area of dust in space, and applying the Mie light-scattering theory to estimate 
the temperatures and emissivities of dust particles in broad size and heliocentric distance ranges. In particular, the model of the 
interplanetary dust cloud by Kelsall et al. (1998, ApJ 508: 44-73), the five populations of interplanetary meteoroids of Divine (1993, 
JGR 98(E9): 17,029-17,048) and the Interplanetary Meteoroid Engineering Model (IMEM) by Dikarev et al. (2004, EMP 95: 109- 
122) are used in combination with the optical properties of olivine, carbonaceous and iron spherical particles. The Kelsall model has 
been widely accepted by the Cosmic Microwave Background (CMB) community. We show, however, that it predicts the microwave 
emission from interplanetary dust remarkably different from the results of application of the meteoroid engineering models. We make 
maps and spectra of the microwave emission predicted by the three models assuming variant composition of dust particles. Predictions 
can be used to look for the emission from interplanetary dust in CMB experiments as well as to plan new observations. 
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1. Introduction 


The unprecedented high-precision surveys of the mi¬ 
crowave sky performed recently by the Wilkinson Microwave 
Anisotropy Probe (WM AP, Bennett e t alJ 1201 3l) and Planck 
llPlanck Collaboration Tll20141) observatories in search and char¬ 
acterisation of the large- and small-scale structure of the Cosmic 
Microwave Background (CMB) anisotropies pose new chal¬ 
lenges for simulation and subtraction of the foreground emission 
sources, including the Solar-system dust. Previous^ templates 
designed for this purpose were based on the iKelsall et al.l 
(1998) model. They indicated little significance of the i nter- 
planetary dust for the WMAP survey (iSchlegel et al.lfl998l) . yet 
remarkable contribution was_ detected at high frequencies of 
Planck (iPlanck Collaboration XIVlf2014l) . 

The angular power spectrum of CMB anisotropies is in 
good agre ement with the inflationary A-cold dark matter model 
([Planck Collaboration B2014I) . However, the reconstructed CMB 
maps at the largest angular scales reveal some intriguing anoma¬ 
lies, among them unexpected alignments of multipole mo¬ 
ments, in particular with directions singled out by the So¬ 
lar system llSchwarz et al.ll20Q4l) . The quadrupole and octopole 
are found to be mutually aligned and they define axes that 
are unusually perpendicular to the ecliptic pole and paral¬ 
lel to the direction of our motion with respect to the rest 
frame of the CMB (the dipole direction). For reviews see 


iBennett et al. ([201 1 [), Coni et al. (1 201 Ol) a ndfeennett et al. ([20 1 3l) 
and Planck Co llaboratio n XXIliP i 20141) as well as Coni et akl 
d201 3T for a detailed analysis of final WMAP and first Planck 
data. It has been suggested that an unaccounted observation bias 
may persist, e.g. yet another foreground source of the microwave 


emission bound to the Solar system. [Dikarev et al.l ( 2009 ) ex¬ 
plored the possibility that the Solar-system dust emits more in¬ 
deed than it was previously thought, and found that the macro¬ 
scopic (> 0.1 mm in size) meteoroids can well contribute ~ 
10 pK in the microwaves, i.e. comparable with the power of the 
CMB anomalies, without being detected in the infrared (IR) light 
and in c luded in the I R-bas ed m odels like tha t of Kelsall et al.l 
(119981) . iBabich et al.l (120071) and lHansen et al.l (l2012h have also 
studied possible contribution of the Kuiper belt objects to the mi¬ 
crowave foreground radiation. [Dikarev et al.l (feOOSl) constructed 
and tested against the WMAP data some dust emission templates. 


Here we improve and extend preliminary estimates of 
[Dikarev et al .i! (12009b . In addition to the Kelsall model, we use 
two meteoroid engineering models to make accurate and thor¬ 
ough maps of the thermal emission from the Zodiacal cloud in 
the wavelength range from 30 pm to 30 mm. 

The paper is organized as follows. Section[2]introduces three 
models of the interplanetary meteoroid environment that we use 
to predict the thermal emission from the zodiacal dust cloud. A 
detailed description of the theory and observations incorporated 
in each model is provided, it may be helpful not only for cosmol- 
ogists interested in understanding the solar-system microwave 
foregrounds, but also for developers of new meteoroid models 
willing to comprehend earlier designs. Section [3] describes the 
thermal emission models. Empirical models as well as the Mie 
light-scattering theory are used to calculate the temperatures and 
emission intensities of spherical dust particles composed of sil¬ 
icate, carbonaceous and iron materials, for broad ranges of size 
and heliocentric distance. The models of the spatial distribution 
of dust and thermal emission are combined in Section[4]in order 


Article number, page 1 oil 171 
































































A&A proofs: manuscript no. paper-con-aa-2 


Table 1 . Frequencies v, wavelengths A and bandwidths Av/v of 
the CIB/CMB observations with COBE DIRBE wavelength bands 4 
through 10, WMAP, and Planck. 


Instrument 

V 

A 

Av/v 

COBE 

61.2 THz 

4.9 pm 

0.13 

DIRBE 

25.0 THz 

12 pm 

0.54 


12.0 THz 

25 pm 

0.34 


5.00 THz 

60 pm 

0.46 


3.00 THz 

100 pm 

0.32 


2.14 THz 

140 pm 

0.28 


1.25 THz 

240 pm 

0.40 

Planck 

857 GHz 

350 pm 

0.25 

HFI 

545 GHz 

550 pm 

0.25 


353 GHz 

850 pm 

0.25 


217 GHz 

1.4 mm 

0.25 


143 GHz 

2.1 mm 

0.25 


100 GHz 

3.0 mm 

0.25 

Planck 

70 GHz 

4.3 mm 

0.2 

LFI 

44 GHz 

6.8 mm 

0.2 


30 GHz 

10.0 mm 

0.2 

WMAP W 

94 GHz 

3.2 mm 

0.22 

V 

61 GHz 

4.9 mm 

0.23 

Q 

41 GHz 

7.3 mm 

0.20 

Ka 

33 GHz 

9.1 mm 

0.21 

K 

22 GHz 

13.6 mm 

0.24 


to make all-sky maps and spectra of the thermal emission from 
the interplanetary dust. Conclusions are made in Sect. 0 

Throughout this paper, we will be dealing with the wave¬ 
lengths mostly, since the dust optical properties are naturally de¬ 
scribed in terms of linear scales. The CMB community, however, 
is more accustomed to frequencies. A conversion table is there¬ 
fore useful for the observation wavebands of infrared detectors 
and radiometers (Table [!}• 


2. Meteoroid Models 

In this paper, we use three recen t and contemporary models 
of the Zodiacal dust cloud due to iKelsall et al.l (11998k . iDivinel 
(fl993h and lDikarev et alJ (120041) . The first model is constrained 
by the infrared observations only from an Earth-orbiting satel¬ 
lite COBE. For brevity, it is referred to as the Kelsall model 
hereafter. The Kelsall model is most familiar to and most of¬ 
ten used by the Cosmic Microwave Background (CMB) research 
community as it helps, and was designed to assess and mitigate 
the contamination of the background radiation maps by the fore¬ 
ground radiation from interplanetary dust. The two other mod¬ 
els incorporate data obtained by diverse measurement techniques 
and serve to predict meteoroid fluxes on spacecraft and to assess 
impact hazards. Being applied in spacecraft design and analysis 
mostly, they are often dubbed as meteoroid engineering models. 
We choose the Divine model and IMEM, with the latter abbre¬ 
viation standing for the Interplanetary Meteoroid Engineering 
Model. We also check if the most recent NASA meteoroid en¬ 
gineering model, MEM (iMcNamara et al.l[2004l) . can be used in 
our study, and explain why it cannot be. 

The Kelsall model does not consider the size distribution of 
particles in the Zodiacal cloud. The size distribution is implic¬ 
itly presented by the integral optical properties of the cloud. In 
contrast, the meteoroid engineering models are obliged to spec¬ 
ify the fluxes of meteoroids for various mass, impulse, or other 


level-of-damage thresholds, hence they provide the size distribu¬ 
tion explicitly. 

Profiles of particle cross-section area density per unit volume 
of space are plotted for the Kelsall and Divine models as well 
as IMEM in Fig. Q] Interestingly, the Kelsall model is sparser 
than both the Divine model and IMEM at most heliocentric dis¬ 
tances. The meteoroid engineering models predict substantially 
more dust than Kelsall does, especially beyond 1 AU from the 
Sun: Divine’s density is almost flat in the ecliptic plane between 
1 and 2 AU, exceeding Kelsall’s density by a factor of 3 in the as¬ 
teroid belt. The IMEM density remains similar to that of the Kel¬ 
sall model up to ~ 3 AU, showing a local maximum beyond the 
asteroid belt near 5 AU from the Sun, in the vicinity of Jupiter’s 
orbit. IMEM possesses the latitudinal distribution somewhat nar¬ 
rower than that of the Kelsall and Divine models. 

All these distinctions stem from different observations and 
physics incorporated in different models. When describing them 
in subsequent subsections, we do not aim at selecting the best 
model. We take all three of them instead, with the intention to 
“bracket” the more complex reality by three different perspec¬ 
tives from different “points of view”. 

2.1. The Kelsall model 

A concise analytical description of the infrared emission from 
the interplanetary dust captured by the COBE Diffuse InfraRed 
Background Experiment (DIRBE) was proposed by Kelsall et al.l 
(11998b . It recognizes five emission components, a broad and 
bright smooth cloud, three finer and dimmer dust bands, and 
circumsolar dust ring along the Earth orbit with an embedded 
Earth-trailing blob. Each component is described by a parame- 
terised empirical three-dimensional density function specifying 
the total cross-section area of the component’s dust particles per 
unit volume of space (Fig. [2}, and by their collective light scatter¬ 
ing and emission properties such as albedo, absorption efficien¬ 
cies, etc. The particle size distributions of the components are 
not considered. Most of the light scattering and emission prop¬ 
erties are neither adopted from laboratory studies of natural ma¬ 
terials nor predicted by theories of light scattering, they are free 
parameters of the model fit to the DIRBE observations instead. 

The smooth cloud is the primary component of the Kelsall 
model. Its density function is traditionally separated into radial 
and vertical terms. The radial term is a power law 1 //?" with 
a = 1.34 ± 0.022 being the slope of the density decay with 
distance from the cloud’s centre R c . The slope is known to be 
equal to one for the dust particles in circular orbits migrating to¬ 
ward the Sun under the Poynting-Robertson drag (orbiting the 
Sun, the particles absorb its light coming from a slightly for¬ 
ward direction due to aberration, and the radiation pressure force 
has a non-zero projection against the direction of their motion; 
that projection causes gradual loss of the orbital energy and de¬ 
crease of the semimajor axis of particle’s orbit). The slope is 
greater than one if the particle orbits are initially eccentric, or if 
the zodiacal cloud is fed from sources broadly distributed over 
radial distances, so that the inner circles of the cloud are sup ¬ 
plied from more dust source s than the outer circles dLeinert et al.l 
11 983HGor’kavvi et al III 997h . 

The vertical term of smooth cloud’s density is represented 
by a widened, modified fan model. The smooth cloud’s symme¬ 
try plane is inclined off the ecliptic plane. This happens because 
the Earth has no influence on the orbital dynamics of the vast 
majority of cloud’s particles. The giant planets, Jupiter primar¬ 
ily, perturb the orbits of dust particles as well as their sources, 
or parent bodies, such as comets and asteroids, and control the 
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Fig. 1. Profiles of particle cross-section area density predicted with selected meteoroid models and compared with each other: ecliptic radial (left) 
and latitudinal at 1 AU from the Sun (right, the vernal equinox is at latitude 0°). 


inclinations of the orbits of the particles. An offset of the center 
of the cloud from the Sun is also allowed in the Kelsall model 
(of the order of 0.01 AU). The vertical motion of the Earth with 
respect to the cloud’s symmetry plane leads to the modulations 
of the infrared emission from the zodiacal cloud reaching 30% 
in the ecliptic polar regions. 

The dust bands are the remnants of collisional disruptions of 
asteroids that resulted in formation of the families of asteroids 
in similar orbits and remarkable dust belts. The dust bands were 
discovered first on the high-resolution maps made by the IRAS 
satellite dLow et akl 1 19841) and re c ogni zed later in the COBE 
DIRBE data as well. iKelsall et alJ (1 19981) introduced three dust 
bands in their model, attributed to the Themis and Koronis aster¬ 
oid families (±1?4 ecliptic latitude), the Eos family (±10°), and 
the Maria/Io family (±15°). The bands have symmetry planes 
different from that of the smooth background cloud. The bands 
are all double, with the density peaking above and below the 
respective symmetry plane, since their constituent particles re¬ 
tain the inclination of the ancestor asteroid orbit but the longi¬ 
tudes of nodes get randomized by the planetary perturbations. 
The vertical motion of such particles is similar to that of an os¬ 
cillator, which moves slower and spends more time near its ex¬ 
tremities, thus the density at the latitudes of each family’s or¬ 
bital inclination is high. One would expect a similarly shaped 
“edge-brightened” radial density with the peaks at the perihe¬ 
lion and aphelion distances of the parent body orbit (e.g. Sect. 3 
in iGor’kavvi et al.l Il997h . However, the dust in asteroid bands 
is proven to be migrati ng toward the Sun under the Poynting- 
Robertson drag: IReachl (119921 ) has demonstrated that the tem¬ 
peratures of particles and parallaxes of the bands measured from 
the moving Earth are both higher than those supposed to be in 
the respective asteroid families of their origin, implying that the 
bands extend farther t owards the Sun than their ancestor bodies. 
Unlike IReachl (119921) . IKelsall et al.1 dl998l) allowed for a partial 
migration only of dust, by introducing a cut-off at a minimal 
heliocentric distance either defined or inferred individually for 
each band. They also fixed rather than inferred most of the dust 
band shape parameters. It is noteworthy that lKelsall et alJ(l!998l) 
did not introduce an explicit outer cut-off for the densities of the 
smooth cloud and dust bands. Instead, they integrated the densi¬ 
ties along the lines of sight from the Earth to the maximal helio¬ 
centric distance of 5.2 AU (close to the Jovian orbital radius). 


The dust particles on nearly circular orbits migrating under 
the Poyinting-Robertson drag toward the Sun are temporarily 
trapped in a_mean-motion resonance with the Earth near 1 AU 
(iDermott et al.lfl984}) . The resulting enhancement of the zodia¬ 
cal cloud density is described in the Kelsall model by the solar 
ring and Earth-trailing blob. The trailing blob is the only com¬ 
ponent of the model that is not static: it is orbiting the Sun along 
a circular orbit with a period of one year, and as the name sug¬ 
gests, its density peak is located behind the moving Earth. The 
solar ring and trailing blob have density peak distances slightly 
outside the Earth orbit, their symmetry planes are inclined off the 
ecliptic plane, and the trailing blob orbits the Sun at a constant 
velocity whereas the Earth velocity is variable due to the eccen¬ 
tricity of our planet’s orbit. Consequently, the Earth moves with 
respect to the solar ring and trailing blob over the course of the 
year much like the other components of the Kelsall model. 

The DIRBE instrument performed a simultaneous survey of 
the sky in 10 wavelength bands at 1.25, 2.2, 3.5, 4.9, 12, 25, 
60, 100, 140, and 240 /-im. It was permitted to take measure¬ 
ments anywhere between the solar elongations of 64° and 124°. 
As the Earth-bound COBE observatory progressed around the 
Sun, DIRBE took samples of the infrared background radiation 
from all over the sky. The foreground radiation due to interplane¬ 
tary dust was a variable ingredient of the sky brightness because 
it depends on the changing position of observatory with respect 
to the cloud. Thus the fitting technique was based on minimizing 
the difference between the brightness variations in time observed 
along independent lines of sight and those predicted by the Kel¬ 
sall model for the same lines of sight, ignoring the underlying 
photometric baselines contaminated or even dominated by the 
background sources. 

The accuracy of the Kelsall model in describing the infrared 
emission from interplanetary meteoroids in DIRBE’s wave¬ 
length bands and within the range of permitted solar elonga¬ 
tions is reported to be better than 2%. Interpolations of the model 
brightnesses between the instrument wavelengths and extrapola¬ 
tions beyond the wavelength and solar elongation ranges may be 
prone to higher uncertainties. Indeed, the concise description of 
the zodiacal emission contribution to the DIRBE data provides 
no clue as to how the light scattering and emission properties of 
its components depend on wavelength between and beyond the 
ten DIRBE wavelength bands. 
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Fig. 2. Radial and latitudinal profiles of particle cross-section area density in the Kelsall model, componentwise. 


When applying the Kelsall model in the far infrared wave¬ 
lengths and microwaves, one should bear in mind that a I re tidy 
iKelsall et al.l (119981) found the 140 and 240 pm bands nearly 
useless in constraining the density distribution parameters due 
to relatively high detector noise and small contribution of the ra¬ 
diation from the interplanetary dust at these wavelengths. Con¬ 
sequently, the inferred parameters of the density distributions in 
the Kelsall model are based on the observations at the wave¬ 
lengths up to 100 pm. The dust particles are efficient in interact¬ 
ing with the electromagnetic radiation if their sizes are not too 
small with respect to the wavelength (2ns > A, where s is the 
particle radius and A is the wavelength). Dust grains with radii 
fro m ~ 1 6 pm are therefore visible at the wavelength of 100 pm. 
iDikarev et al .l (12009b argue that a considerable amount of mete¬ 
oroids are present in the Solar system with the sizes bigger than 
~ 100 pm which are outshined in the infrared light by more 
abundant and ubiquitous smaller dust grains. The lo ngest obser¬ 
vation wavelength of DIRBE at which lKelsall et al.l l 1 998 ) could 
still constrain the density distributions of dust is short of being 
capable to resolve these particles unless they compose distinctive 
features like dust bands. 

We use the formulae for the density functions provided by 
Planck Collaboration XlM dm! since the original paper by 
Kelsall et al.l 1199a ) contains typographical misprints in the def¬ 
initions of the asteroid dust band and circumsolar ring densities. 

2.2. The Divine model 

A model of the interplanetary meteoroid environment ( Divine! 
119931) constructed to predict fluxes on spacecraft anywhere in 
the Solar system from 0.05 to 40 AU from the Sun, is composed 
of five distinct populations, each of which has mathematically 
separable distributions in particle mass and in orbital inclination, 
perihelion distance, and eccentricity (Fig. |3j. Using the distribu¬ 
tions in orbital elements. Divine explicitly incorporated Keple- 
rian dynamics of meteoroids in heliocentric orbits. 

The Divine model is constrained by a large number of diverse 
meteoroid data sets. The mass distribution of meteoroids from 
1CU 18 to 10 2 g is fitted to the interplanetary meteoroid flux model 
dGriin et alii 19851) . which in turn is based on the micro-crater 
counts on lunar rock samples retrieved by the Apollo missions, 
i.e. the natural surfaces exposed to the meteoroid flux near Earth, 


and data from several spacecraft. The orbital distributions are de¬ 
termined using meteor radar data, observations of zodiacal light 
from the Earth as well as from the Helios and Pioneer 10 inter¬ 
planetary probes, and meteoroid fluxes measured in-situ by im¬ 
pact detectors on board Pioneer 10 and 11, Helios 1, Galileo and 
Ulysses spacecraft at the heliocentric distances ranging from 0.3 
to 18 AU. The logarithms of the model-to-data ratios were mini¬ 
mized in a root-mean-square sense. When modeling the zodiacal 
light. Divine assumed that the scattering function is independent 
of meteoroid mass. The albedos of dust particles were defined 
somewhat arbitrarily in order to hide the populations necessary 
to fit the in-situ measurements from the zodiacal light observa¬ 
tions. No infrared observations were used to constrain the model. 

The core population is the backbone of the Divine model, it 
fits as much data as possible with a single set of distributions. 
The distribution in perihelion distance r n of the core-population 
particles can be used as strictly proportional to r^ 13 up to 2 AU 
from the Sun. This function resembles a spatial concentration, 
and its slope is close to Kelsall’s a = 1.34 for the radial density 
term: much like the infrared data from COBE , the zodiacal light 
observations from the Earth and interplanetary probes demand 
that the slope be steeper than a unit exponent. The eccentricities 
are moderate (peaking at 0.1 and mostly below 0.4), inclinations 
are small (mostly below 20°). 

The other populations fill in the gaps where one population 
with separable distributions is not sufficient to reproduce the ob¬ 
servations. Their names hint at distinctive features of their orbital 
distributions. The inclined population is characterized by incli¬ 
nations largely in the range 20°-60°, its eccentricities are all be¬ 
low 0.15. In contrast, the eccentric population is composed of 
particles in highly eccentric orbits (the eccentricities are largely 
between 0.8 and 0.9), the inclinations are same as in the core 
population. These two populations are added in order to com¬ 
pensate for a deficit of meteoroids from the core population with 
respect to in-situ flux measurements on board the Helios space¬ 
craft. (The eccentric population is also used to match the inter¬ 
planetary flux model between 1CT 18 and 10 -15 g.) The halo popu¬ 
lation has a uniform distribution of orbital plane orientations and 
“surrounds” the Sun as a halo between roughly 2 and 20 AU. It 
patches the meteoroid model where the above-listed populations 
are not sufficient to reproduce the Ulysses and Pioneer in-situ 
flux measurements. The eccentric, inclined and halo populations 
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Fig. 3. Radial and latitudinal profiles of particle cross-section area density in the Divine model, componentwise. 


are composed of meteoroids much smaller than those significant 
enough for predictions of the infrared and microwave observa¬ 
tions. 

The asteroidal population is described last but it makes a 
huge difference between Kelsall’s and Divine’s models. Already 
when fitting the first population to the meteor radar data, Divine 
observed that the results for the distribution in perihelion dis¬ 
tance show a minimum near r„ = 0.6 AU, suggesting that two 
components are involved (i.e., inner and outer). The radial slope 
of the inner component was consistent with that demanded by 
the zodiacal light observations, and the inner and outer compo¬ 
nents contributed 45% and 55%, respectively, to the flux near 
Earth. The outer component’s concentration peaked outside the 
Earth orbit, in the asteroid belt. Since the main contribution to 
the zodiacal light comes from the particles with masses smaller 
than 1CT 5 g and the median meteoroid mass of the meteor radar 
was 1CT 4 g, it was reasonable to ascribe the inner fraction to a 
population dominated by smaller particles responsible for the zo¬ 
diacal light, and the outer fraction to another one dominated by 
larger particles detected as meteors only with the radar. The two 
fractions were divided into the core and asteroidal populations, 
accordingly. Their inclination and eccentricity distributions are 
identical, but the distributions in perihelion distance and mass 
are different. The asteroidal population is used to fit the inter¬ 
planetary flux model at large masses > 1CT 4 g. 

Derived from data analysis rather than postulated theoreti¬ 
cally, the distinction between the core and asteroidal populations 
in the Divine model separates - both on the mass and perihe¬ 
lion distance scales - the small dust particles migrating toward 
the Sun under the Poynting-Robertson drag, the concentration of 
which increases with the decrease of the heliocentric distance, 
and their immediate parent bodies, the larger particles swarming 
further away from the Sun, i.e. in the asteroid belt according to 
iDivind (ll 9931) . 

2.3. I MEM 

The Interplanetary Meteoroid Engineering Model (IMEM) is 
developed for ESA by [Dikarev et al. ( 12004 !). Like the Divine 
model, it uses the distributions in orbital elements and mass 
rather than the spatial density functions of the Kelsall model, en¬ 
suring that the dust densities and fluxes are predicted in accord 


with Keplerian dynamics of the constituent particles in heliocen¬ 
tric orbits. Expansion of computer memory has also enabled the 
authors of IMEM to work with large arrays representing multi¬ 
dimensional distributions in orbital elements to replace Divine’s 
separable distributions. The distribution in mass remains separa¬ 
ble from the 3-D distribution in inclination, perihelion distance 
and eccentricity. 

IMEM is constrained by the micro-crater size statistics col¬ 
lected from the lunar rocks retrieved by the Apollo crews, COBE 
DIRBE observations of the infrared emission from the inter¬ 
planetary dust at 4.9, 12, 25, 60, and 100 pm wavelengths, and 
Galileo and Ulysses in-situ flux measurements. An attempt to in¬ 
corporate new meteor radar data from the Advan ced Meteor Or¬ 
bit Radar {AMOR, see Galligan & Baggaievll2004l) was not suc¬ 
cessful jDikarev et al.l 120051) . as the latitudinal number density 
profile of meteoroids derived from the meteor radar data stood 
in stark contrast to that admissible by the COBE DIRBE data, 
due to incomplete understanding of the observation biases. 

The meteoroids in IMEM are distributed in orbital ele¬ 
ments in accord with an approximate physical model of the ori¬ 
gin and orbital evolution of particulate matter under the plan¬ 
etary gravity, Poynting -Robertson effect and mutual collisions. 
[Grim & Dikarev! ( 2009 ) visualise the IMEM distributions of me¬ 
teoroids in mass and in orbital elements in graph and color plots. 
Cross-section area density profiles are shown here in Fig. [4] 

All interplanetary meteoroids are divided into populations 
by origin/source and mass/dynamical regime. IGriin et all ( 1985 ) 
constructed a model of the mass distribution of interplane¬ 
tary meteoroids. Based on this model, they estimated particle 
lifetimes against two destructive processes, mutual collisions 
and migration downward toward the Sun due to the Poynting- 
Robertson drag (terminated by particle evaporation or thermal 
break-up). Figure [5] shows that the lifetime against mutual colli¬ 
sions is shorter than the Poynting-Robertson time for the mete¬ 
oroids bigger than ~ 10 s g in mass, located at 1 AU from the 
Sun. The destruction of meteoroids is dominated by a factor of 
ten in rate by the respective process already one order of magni¬ 
tude away from this transition mass in either direction. The rates 
vary with the distance from the Sun, however, the transition mass 
rema ins nearly unchange d (Fig. 6 in IGriin et al.lll985h . That is 
why [Dikarev et al.l d2004l) could simplify the problem by divid¬ 
ing the mass range into two subranges with distinct mass dis- 
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Fig. 4. Radial and latitudinal profiles of particle cross-section area density in IMEM, componentwise. The density plots are not smooth since the 
orbital distributions in IMEM are approximated by step functions. 


tribution slopes and dynamics. Meteoroids below the transition 
mass of 10~ 5 g are treated as migrating from the sources toward 
the Sun under the Poynting-Robertson drag, meteoroids above 
the transition mass are assumed to retain the orbits of their ul¬ 
timate ancestor bodies over their short lifetimes until collisional 
destruction. 

The ultimate ancestor bodies of interplanetary dust in IMEM 
are comets and asteroids. A catalogue of asteroids is used to 
account for the latter source of meteoroids. Three distinct pop¬ 
ulations of asteroids are recognized and allowed to have inde¬ 
pendent total production rates: the Themis and Koronis families 
(semi-major axes 2.8 < a < 3.25 AU, eccentricities 0 < e < 0.2, 
and inclinations 0 < i < 3?5), Eos and Veritas families (2.95 < 
a < 3.05 AU, 0.05 < e < 0.15, 8?5 < i < 11?5), and the main 
belt (a < 2.8 AU). The orbital distributions of meteoroids more 
than 10~ 5 g in mass are constructed by counting the numbers 
of catalogued asteroids per orbital space bins. The orbital distri¬ 
butions of meteoroids less than 10 s g in mass are described as 
the flow of particles along the Poynting-Robertson evolutionary 
paths dGor’kavvi et al.li 1997b starting from the source distribu¬ 
tions defined earlier. The mass distributions are adopted piece- 
wise from the interplanetary meteoroid flux model bv lGriin et al.1 
d 1985b as shown in Fig. [5] 

The orbital distributions of meteoroids from comets cannot 
be defined as easy as those of meteoroids from asteroids. Be¬ 
cause of a number of loss mechanisms, such as comet nuclei de¬ 
cay, accretion, tidal disruption or ejection from the Solar system 
by planets, very few comets are active at a given epoch and listed 
in the catalogues. Moreover, the catalogues are prone to obser¬ 
vation biases since the comet nuclei are revealed by gas and dust 
shed at higher rates at lower perihelion distances. 

In order to de scribe the orbital distributions of meteoroids of 
cometary origin, iDikarev & Griinl ( 2004 ) proposed an approx¬ 
imate analytical solution of the problem of a steady-state dis¬ 
tribution of particles in orbits with frequent close encounters 
with a planet. This solution is applied to represent the orbital 
distributions of meteoroids from comets in Jupiter-crossing or¬ 
bits, i.e. the vast majority of catalogued comets. Big meteoroids 
with masses greater than ICE 5 g are confined to the region of 
close encounters with Jupiter. There is only one quantity that is 
approximately conserved in the region of close encounters with 


Jupiter. It is the Tisserand parameter T 

Cl] / Cl _ T 

T = — + 2 . /—(1 - e 2 )cos i, 
a V flj 

with aj being the semimajor axis of Jovian orbit. The num¬ 
bers of meteoroids in T -layers are therefore stable over a long 
period of time. They are free parameters to be determined 
from the fit of IMEM to observations. The relative distribu¬ 
tions in orbital elements with in each T -layer are found theoreti¬ 
cally (IDikarev & Griinll2004h . 

Finer dust grains with masses less than 1CT 5 g leak from 
this region into the inner solar system and then migrate to¬ 
ward the Sun under the Poynting-Robertson drag. The Poynting- 
Robertson migration time is much longer than the close approach 
time in the region of close encounters with Jupiter, consequently, 
the orbital distributions of all meteoroids are shaped there by 
the gravitational scattering on Jupiter. It is only the dust leak¬ 
ing through the inner boundary of the region that is distributed 
as the flow of particles along the Poynting-Robertson evolution¬ 
ary paths (iGor’kavvi et al.lfl997t) . Their distributions are also es¬ 
tablished theoretically in IMEM, with the normalisation factors 
being inferred from the fit to observations. 

An empirical finding by Divine of a helpful separation of the 
bulk of meteoroid cloud into the core and asteroidal populations 
- segregation of big and small dust particles - has become one 
of the physical assumptions in IMEM. 

The weights of populations of dust from comets and aster¬ 
oids in IMEM were fitted to in-situ flux measurements using the 
Poisson maximum-likelihood estimator and to infrared observa¬ 
tions using the Gaussian maxiumum-likelihood estimator. 

IMEM was tested against those data not incorporated in 
the model by IDikarev et al.l ( 2005 ). confirming that the or¬ 
bital evolution approach allows for more reliable extrapola¬ 
tions of the observations and measurements incorporated in the 
model, it is compared w ith several other meteoroid models by 
iDrolshagen et ali (120081) . 

2.4. MEM 

The Meteoroid Engin eering Model (MEM) described by 
[McNamara et ali ( 2004 1 is another recent development of the or- 
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Fi fe 5 - Cu mulati ve mass distribution of interplanetary meteoroids 
dCirun et al.l 11985L top), meteoroid lifetimes against mutual collisions 
and Poynting-Robertson drag (middle), and cumulative mass distribu¬ 
tions of meteoroid populations in IMEM dDikarev et aljl2005l . bottom). 
The flux of meteoroids with m > ICE 5 g is lower in IMEM than in the 
Grim model since the impact velocities of big meteoroids are higher at 
1 AU from the Sun in IMEM and they produce larger craters on the lunar 
rocks than in the Griin model, which assumed a single impact velocity 
for all meteoroid sizes. 


bital distributions and software to predict fluxes on spacecraft in 
the Solar system and near Earth. It is constrained by the Earth- 
based meteor radar data CMOR (Canadian Meteor Orbit Radar) 
and zodiacal light observations from two interplanetary probes, 
Helios 1 and 2. The nominal heliocentric distance range at which 
the model is applicable, from 0.2 to 2 AU, is rather limited, how¬ 
ever. Even though the meteoroid mass range of 10 6 to 10 g cov¬ 
ered by the model is extremely interesting for the purposes of 
our study, some assumptions made by the authors are rather ar¬ 
guable. In particular, the mutual collisions between meteoroids 
are considered in the model ignoring dependence of the mete¬ 
oroid collision probability on its size, w h e re as iGriln et aTl ( 1985 ) 
calculated that the lifetime against collisional disruption of me¬ 
teoroids varies by a factor of 100 between the masses of 10 6 and 
1 g (their Fig. 6 and our Fig. [5] in the text above)! It is not only 
the mass distribution of particles that is determined by the colli¬ 
sion probability, but also the orbital distributions: as the lifetime 


against collisions becomes shorter than the Poynting-Robertson 
time by orders of magnitude, the meteoroids can no longer sur¬ 
vive a travel from their sources toward the Sun and the orbital 
distributions of bigger particles are more similar with the dis¬ 
tributions of their sources than the distributions of smaller par¬ 
ticles. MEM is unfortunately unable to reveal and describe this 
important feature of the interplanetary meteoroid cloud. 

2.5. The Cross-Section Area Distributions 

Figure [6] maps the total cross-section area of dust particles per 
unit volume of space, as a function of the position in the Solar 
system, for each component or population of every meteoroid 
model introduced in the previous Section. The X and Z coordi¬ 
nates are measured from the Sun, with the X axis pointing to 
the vernal equinox and Z axis being perpendicular to the eclip¬ 
tic plane. The integral of the cross-section area density along a 
line of sight gives the optical depth of the cloud along that line 
of sight, assuming the geometric-optics approximation, spherical 
particles, and ignoring that the particle efficiencies in absorbing 
and scattering light can in fact be higher or lower than unity. 

The Kelsall model is reproduced in the upper row of maps. 
Three dust bands are shown together in the middle plot. The 
smooth cloud is to the left, and the solar ring is to the right of 
it. The dominant component of the Kelsall model, the smooth 
cloud is getting higher in density only toward the Sun and the 
symmetry plane slightly inclined off the ecliptic plane. The dust 
bands are the only component of the model bulking beyond the 
Earth orbit, and they are by definition bound to the ecliptic lati¬ 
tudes of their ancestor asteroid families’ orbit inclinations. Both 
the smooth cloud and the dust bands extend up to the heliocen¬ 
tric distance of 5.2 AU only, since iKelsall et al. ( 1998 ) stopped 
integrating the model densities there. A more rigorous model of 
the bands composed of dust migrating toward the Sun due to 
the Poynting-Robertson effect would put their densities’ cut-offs 
at the outer boundaries of the corresponding asteroid families. 
Note, however, that most of the thermal infrared emission from 
the interplanetary dust observed by COBE was due to the dust 
within 0.5-1 AU from the Earth orbit, since COBE could not be 
pointed too close to the Sun and dust is too cold and inefficient 
at thermal emission far from the Sun dDikarev et al.ll2009L their 
Sect. 2.1 and Fig. 2). Thus for the purpose of modelling the in¬ 
frared thermal emission observed from the Earth, the behaviour 
of the density at multiples of an astronomical unit was not im¬ 
portant. This may not be the case for the microwave emission 
though. 

The next two rows of maps depict the populations of the Di¬ 
vine model. The core and asteroidal populations are important 
in the infrared and microwave emission ranges. They are dis¬ 
tributed remarkably different in space, with the core population 
density growing higher toward the Sun and the asteroidal popula¬ 
tion density peaking beyond the orbit of the Earth. The two popu¬ 
lations are also composed of particles of different sizes: the core 
population mostly smaller and the asteroidal population mosly 
bigger than ~ 50 /mi. The halo, inclined and eccentric popula¬ 
tions are provided for the sake of completeness. Their tiny, up to 
~ 10 /mi-sized particles are ignorable in the wavelength ranges 
of interest. 

The bottom row of maps and the right map in the third row 
exhibit the populations from IMEM. The asteroidal dust parti¬ 
cles bigger than 1CU 5 g in mass are confined to the asteroid belt, 
naturally. They are not migrating from their ancestor body orbits 
toward the Sun nor extend beyond the outer edge of the asteroid 
belt. The smaller dust particles move their ways toward the Sun 
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due to the Poyting-Robertson effect. The inclinations of their or¬ 
bits are intact, since IMEM does not take into account any plan¬ 
etary perturbations other than gravitational scattering by Jupiter 
that require a close approach to the giant planet. Consequently, 
the latitudinal distribution of their density is independent on he¬ 
liocentric distance away from the source region, i.e. asteroid belt. 

Big cometary particles have a density peak along the orbit 
of Jupiter. This occurs because all their orbits are required to 


approach Jovian orbit within 0.5 AU or less. Small cometary 
particles leak through the inner boundary of the region of close 
encounters with Jupiter and migrate toward the Sun under the ac¬ 
tion of the Poynting-Robertson effect, they reach the highest den¬ 
sity at the shortest heliocentric distances. A cometary dust den¬ 
sity enhancement in the form of a spherical shell with a radius of 
5.2 AU is a negligibly small defect caused by the assumption of 
a uniform distribution of particle orbits in longitudes of perihe- 
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lia in IMEM. Concentric spherical “shells” of different densities 
on the map of cometary meteoroids with m > 10 5 g are due to 
finite discretization of the distribution in perihelion distance. 

Cumulative distribution functions of meteoroid cross-section 
area in mass and size for e ach population of the m eteoroid mod¬ 
els by IDivinel (1 19931) and iDikarev et alJ (|2004) are plotted in 
Fig. □ (normalized to unity). Note that one population of the 
Divine model has an assumed material density of 0.25 g cm -3 , 
whereas the standard value assumed for all other populations as 
well as in many other meteoroid models is 2.5 g cm~ 3 . 

The Kelsall model bets on essentially single-population rep- 
resentaion of the interplanetary meteoroid cloud, with the solar 
ring and dust bands being rather minor features (Fig. [2j. The 
meteoroid engineering models, however, state it clearly that the 
cloud is composed of different sorts of dust in major populations 
that are also distributed differently in space. We will see below 
how this affects the microwave emission predictions. 


3. Thermal Emission Models 

The thermal emission from dust particles is expressed with re¬ 
spect to the blackbody emission B\(T) at the wavelength A and 
temperature T, using an emissivity modification factor which 
matches t he absorption efficiency f actor 2abs(s, A) for the parti¬ 
cle size s dBohren & Huffmanlll983l) . gabs is defined as the ratio 
of the effective absorption cross-section area of the particle to its 
geometric cross-section area. 

3 .1. The Kelsall model 

The light scattering and emission properties of dust are not de¬ 
fined in the Kelsall model for the wavelengths between and be¬ 
yond those of the COBE DIRBE instrument. Figure [8] shows 
Kelsall’s emissivity modification factors, or absorption efficien¬ 
cies to preserve the uniformity of terms, for the DIRBE wave¬ 
length bands. iPlanck Collaboration XIVl (120141) used the cross- 
section area density of the Kelsall model and found the absorp¬ 
tion efficiencies for its components to describe approximately 
the thermal emission from interplanetary dust at the wavelengths 
of Planck’s High Frequency Instrument (HFI). iFixsen & Dwekl 


d2002l) used the FIRAS (Far-Infrared Absolute Spectrometer) 
data from COBE and found that the annually averaged spec¬ 
trum of the zodiacal cloud can be fitted with a single blackbody 
at a temperature of 240 K with an absorption efficiency being 
flat at wavelengths shorter than 150 jum and gabs K 4 2 beyond 
150 fim. 


The absor ption efficiencies of dust from all three bands co¬ 
incide in Kelsall et al. (19981) for the wavelengths up to 240 //m. 
[Planck Collaboration XIVl d2014l) have removed this constraint 
and allowed for individual weights of the band contributions in 
the microwaves. They found that the bands keep high <2abs ~ 1 
up to A ~ 3 mm, implying that their constituent particles are 
macroscopic. The absorption efficiency of the smooth cloud de¬ 
cays in the micro waves, however, not exa ctly as sharp as a simple 
approximation of lFixsen & Dwekj ( 2002 ) suggests. 


The inverse problem_solution sometimes led 

IPlanck Collaboration XIVl d2014l) to negative absorption 
efficiencies of the smooth cloud or circumsolar dust ring and 
Earth-trailing blob. (Those negative efficiencies are simply 
missing from Fig. [8] at certain wavelengths, as the logarithmic 
scale does not permit negative ordinates.) Obviously, some 
components of the Kelsall model were used by the fitting 
procedure to compensate for excessive contributions from the 
other components in such cases. This is a strong indication of 
insufficiency of the Kelsall model in the microwaves. 


Our study requires the absorption factors even further in 
the mi crowa ves. We use the numbers i n ferred by iKelsall et al.l 
(1 19981) and [Planck Collaboration XIVl (120141) whenever pos¬ 
sible, i.e. when they are available and positive. A nega¬ 


tive absorption efficiency found for the smooth cloud by 
IPlanck Collaboration XIVl (120141) at A = 2.1 mm is replaced 
with the result of interpolation of the positive efficiencies found 
at A = 1.4 and 3 mm, whereas negative efficiencies for the cir¬ 
cumsolar ring at A = 0.85, 1.4 and 2.1 mm are simply nullified. 
The absorption efficiency is extrapolated beyond the Planck/HFl 
wavelengths (A > 3 mm, i.e. in the WMAP range) using the 
approximation due to iFixsen & Dwekl (f2 0~02li for the smooth 
cloud (2abs K 4 2 . similar to Maris et alJ (120061) who assessed 
the level of contamination of the Planck data by the Zodiacal 
microwave emission based on the Kelsall model, and using flat 
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2abs = 1, 0.5, 1, and 0.1 for the asteroid dust bands 1,2, 3, and 
the circumsolar ring with the trailing blob, respectively. 

3.2. Selected Materials and Absorption Efficiencies 

Predictions of the thermal emission from interplanetary dust 
clouds using the meteoroid engineering models require opti¬ 
cal properties of constituting particles. Following IDikarev et all 
(2009). we use the Mie light-scattering theory and optical con¬ 
stants of silicate and amorphous carbonaceous spherical particles 
in order to estimate the absorption of solar radiation and thermal 
emission by meteoroids. IDikarev et al. ( 200 9!) banned iron from 
substances for a hypothetical meteoroid cloud that could be re¬ 
sponsible for anomalous CMB multipoles. In this paper, how¬ 
ever, we take it back into account since iron is known to com¬ 
pose some meteorites as well as asteroid surfaces, it is present 
in the interplanetary dust particles and therefore it must be con¬ 
sidered when discussing the real, not hypothetical Solar-system 
medium. Water and other ices are still ignored as extremely in¬ 
efficient emitters of the microwaves. 

Dealing with the microwave thermal emission from dust with 
an assumed temperature. IDikarev et al.l(l2009l) did not need opti¬ 
cal constants for the visual and near-infrared wavelengths. Here 
we calculate the temperatures of dust in thermal equilibrium and 
the absorption efficiencies are required for the brightest part of 
the Solar spectrum. Figure [9] plots the data for the wavelengths 
between 0.1 and 100 pm. 

3.3. The Dust Temperatures 

Let us now calculate and discuss the temperatures of dust par¬ 
ticles of different sizes at different distances from the Sun for 
the substances selected in the previous section. Figure [TO] plots 
the equilibrium temperatures of spherical homogeneous particles 
composed of silicate, carbonaceous and iron materials as well as 
the dust temperature used in the Kelsall model. 

The equilibrium temperature is found from the thermal bal¬ 
ance equation 

/-*oo s-»oo 

ns 2 e abs (s, W) d/l = 4ns 2 Q abs (s, A)B a (T d ) d A, 

Jo Jo 
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( 1 ) 


where s is the radius of a spherical dust grain, A is the wave¬ 
length, To is the incident Solar radiance flux, B iiTn) is the 
blackbody radiance at the dust particle’s temperature To- As the 
left-hand side provides the total energy absorbed from a mono- 
directional incident flux, the right-hand side gives the total en¬ 
ergy emitted, omni-directionally. By denoting the absorption ef¬ 
ficiency averaged over the Solar spectrum with Q a , and the same 
quantity averaged over a blackbody spectrum at temperature T d 
with Q(To), then using the Stefan-Boltzmann law and the So¬ 
lar constant, one can rewrite Eq. © in a more concise form (cf. 
lReachl!T98i : 

t d = 279 k [fi0/fi(r D )] 1/4 () 1/2 ’ (2) 


where R is the distance from the Sun. A perfect black body 
with Q a bs = 1 throughout the spectrum has therefore a temper¬ 
ature of 279 K at 1 AU from the Sun, inversely proportional to 
the square root of the distance. This inverse-square-root trend is 
often closely followed by the real dust particles. 

The Kelsall model uses the temperature given by equation 


T = 286 K 


/I AU 

\~R~ 


0.467 


(3) 


iKelsall et alJ ( 1998 1 emphasize that the dust temperature at 1 AU 
and absorption efficiencies could not be determined indepen¬ 
dently, so that the temperature was found by assuming the 
smooth cloud to be the dominant component with its spectrum in 
the mid-IR being that of a pure blackbody (i.e., unit absorption 
efficiencies at A = 4.9, 12, and 25 pm). 

The equilibrium temperatures of dust particles are shown in 
Fig.llOlfor a broad range of heliocentric distances. Temperatures 
used in the Kelsall model are plotted for comparison as well. 
Temperatures of the micrometer-sized particles are in most ob¬ 
vious disaccord with Kelsall’s model, but they tend to be higher 
than temperatures of the larger particles of their composition as 
well. The reason is simple: their absorption efficiencies are too 
low at the wavelengths of ten and more micrometers, at which 
their larger counterparts emit the energy absorbed from the So¬ 
lar radiation flux, thus they warm to higher temperatures in order 
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Fig. 9. Optical constants and abso rption efficiencies for olivine and iron of Polla ck et alJ d!994h. and for amorpho us carbonaceous grains from the 
'BE’ sample of lZubko et alJ dl996h used in this paper, in addition to the data introduced bv TPikarev et al.l 12009 1). The absorption efficiencies are 
calculated for the sizes of spherical dust grains indicated in the bottom row of plots. 


to shift the maximum of the Planck function to shorter wave¬ 
lengths. Iron is an exception from this rule, all particles made 
of this substance are considerably warmer, being unable to emit 
the radiation thermally at longer wavelengths regardless of size. 
(It does not contradict to the fact that communications in ra¬ 
diowaves are assisted by iron and other metallic antennas, as iron 
is very good at scattering the electromagnetic emission at longer 
wavelengths simultaneously!) 

Amorphous carbonaceous particles (labeled “Zubko BE” in 
our plots) from ~ 10 p m in size are in the best agreement with 
the Kelsall model temperatures, among the real species con¬ 
sidered here. This is a good indication for the Kelsall model, 
since most interplanetary meteoroid s are indeed expected to be 
comp osed of carbonaceous material (Reach 19881 iDikarev et aD 
2009). The temperatures of the silicaceous particles are some¬ 
what lower, as olivine is rather transparent in the visual light that 
it absorbs, remaining opaque in the infrared light that it emits, 
but this difference is quite negligible starting from particle radii 
of 100 yum, especially at the heliocentric distances beyond 1 AU 

(Fig. HD. 

In subsequent sections, however, we assume a single temper¬ 
ature for all dust particles of each substance, regardless of size. 
Our results are therefore less precise for dust particles smaller 
than ~ 100 fim in radii. Fortunately for achievement of this pa¬ 
per’s goal, the microwave emission from large particles is more 
important to predict accurately, since small particles are not effi¬ 
cient emitters at that long wavelengths. 


to the Zodiacal cloud, measured in units of power per unit solid 
angle and unit wavelength, observed from the location r in the 
direction specified by a unit vector p, is given by 



Ba(T d (\t + Zp)|)C abs (d, r + Zp)dZ, 


(4) 


where C abs (/l, r) is the absorption cross-section of dust per unit 
volume of space. In the Kelsall model, C abs is a sum of the prod¬ 
ucts of the total cross-section area of dust by the emissivity mod¬ 
ification factors for the wavelength A, over all model compo¬ 
nents. Meteoroid engineering models require an evaluation of a 
nested integral over the particle size s: 

CM r) = - J' tts 2 QM s ) ^ ds, (5) 

where N(s, r) is the number density of meteoroids with radii 
greater than s at the location r. The minus sign in Eq. © is nec¬ 
essary since dN(s,r)/ds < 0. Figure [IT] shows that the equilib¬ 
rium temperature To depends on particle size as well, however, 
this dependence is important only for the small dust grains with 
s < 100 /rm, which can be safely neglected in the microwaves. 

Note that no color correction factor is applied, so that a dif¬ 
ference of at least several percent is to be expected between 
our intensities and Kelsall’s Zodiacal Atlas for the wavelengths 
where a direct comparison is possible. 


4. The Microwave Thermal Emission from Dust 

We are ready now to make maps of the thermal emission from 
the Zodiacal dust for the broadest range of wavelengths and ob¬ 
servation points. Indeed, the surface brightness of the sky due 


4. 1. All-Sky Maps 

All-sky maps of the thermal emission from the interplanetary 
dust are made in Fig.[12]using the Kelsall model as well as me¬ 
teoroid engineering models assuming that their particles are car¬ 
bonaceous. Five observation wavebands of COBE, Planck and 
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Fig. 12. All- sky maps of the therma l emission from the Zodiacal cloud, according to the Kelsall model and two meteoroid engineering mod¬ 
els jDivinel 1993i : iDikarev et all2004l) populated by the carbonaceous particles, as seen from Earth at the fall equinox time, for selected wavebands 
of the COBE , Planck and WMAP observatories. The reference frame is ecliptic, each map's center points to the antisolar direction (i.e., the vernal 
equinox), a 30°-wide band around the Sun is masked. The grey scale of the upper row of maps is in MJy sterad' 1 , the other rows are in pK of a 
temperature in excess of the CMB. Each map has its own brightness scale. 


WMAP are selected to demonstrate the transformation of the 
brightness distribution from the far infrared light (240 pm) to 
microwaves (13.6 mm). 

At short wavelengths, emission from fine dust grains is dom¬ 
inant, and all models of the Zodiacal cloud show the maximum 
emission near the Sun. The smooth background cloud is re¬ 
sponsible for this brightness distribution in the Kelsall model, 
the core population in the Divine model and the meteoroids 
with masses m < ICE 5 g migrating toward the Sun under the 
Poynting-Robertson drag in IMEM are also concentrated near 


the Sun. As the wavelength grows, the fine dust grains fade out, 
but bigger meteoroids remain bright. The components of the Kel¬ 
sall model possessing spectra of macroscopic particles, i.e. the 
asteroid dust bands. Divine’s asteroidal population and the me¬ 
teoroids with masses m > 10~ 5 g in IMEM are located mostly 
beyond the orbit of the Earth. Instead of the brightness peak of 
the fine dust at small solar elongations, the big meteoroids ap¬ 
pear as a band, or bands, along the ecliptic, broadening in the 
antisolar direction: the surface brightness of diffuse objects does 
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Fig. 13. Brightness profiles of the thermal emission from the Zodiacal cloud along the great circles in the ecliptic plane (ecliptic latitude / = 0) and 
perpendicular to the ecliptic plane at solar elongation e = 90°, according to the Kelsall model and two meteoroid engineering models populated by 
the carbonaceous particles, as seen from Earth, for selected wavebands of the COBE, Planck and WMAP observatories. The excess temperatures 
are calculated w.r.t. the CMB. 


not depend on the distance from the observer, but their angular 
size increases due to their proximity in the antisolar direction. 

Figurefl3lshows the brightness profiles built with the models 
from Fig. [12] for Earth-bound observatories scanning the celes¬ 
tial sphere in two important great circles, one in the ecliptic plane 
(zero ecliptic latitude) starting from the vernal equinox counter¬ 
clockwise, and another perpendicular to the ecliptic plane, at the 
solar elongation of 90°, from the ecliptic north pole toward the 
apex of Earth’s motion about the Sun. The scan in the ecliptic 
plane is only unmasked, however, when the target area on the 
sky was visible with WMAP (solar elongations from 90° to 135°) 


or COBE DIRBE (solar elongations from 64° to 124°). Note that 
both profiles match at scan angles of 90° and 270° since the great 
circles are crossing there. 

The Divine model and IMEM with the carbonaceous me¬ 
teoroids appear to be substantially brighter than the Kelsall 
model. Note that IMEM was fitted to the DIRBE data up to 
/l = 100 /rm only, using a composition of dust different from 
the triple used here, and Divine used no infrared data at all. 
Thus no exact match with the Kelsall model is ex pected, espe¬ 
cially at the wavelengths A > 100 jum for which iKelsall et all 
( 1998 ) themselves found the COBE DIRBE data insufficiently 
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Fig. 10. Equilibrium temperatures of homogeneous spherical dust par¬ 
ticles composed of different substances. Fig. 11. Equilibrium temperatures of homogeneous spherical dust par¬ 

ticles composed of different substances. 


constraining the model (Sect. ED> . The difference is substantially 
lower (25-3 0%) at s horter wavelengths (12-60 pm). However, 
IPlanck Collaboration XIVl d2014h inferred the absorption effi¬ 
ciencies of the Kelsall model components for the Planck wave¬ 
bands as well, and reported a lower emission from the Zodiacal 
cloud as shown in Fig. [12] and [13] on the plots for A = 850 pm 
and 2.1 mm. 

To pave the way for an explanation of the discrepancy, we 
remind that the interplanetary dust is not entirely carbonaceous, 
and meteoroids of other chemical composition yield lower in¬ 
frared and microwave thermal emission. Besides, meteoroid en¬ 
gineering models are fitted to many other data sets, typically with 
an accuracy substantially lower than the models of the electro¬ 
magnetic emission from dust. Finally, the fitting of the Kelsall 
model to the Planck data by IPlanck Collaboration XIVl (12014|) 
resulted in confusingly negative absorption efficiencies of dust 
sometimes, implying that the model may be incomplete, with 
some model components used by the fitting procedure to com¬ 
pensate for excessive emission from the other components. One 
can also see that the Kelsall model is deficient in the ecliptic 
plane, with the asteroid dust bands allowed to shine just a lit¬ 
tle above and below that plane only, and not e.g. in the ecliptic 


plane. A new fit to the Planck data using different meteoroid 
models would therefore be of great value, since these models 
may be better in reproducing the relative brightness distribution 
of the microwave emission from dust, but they also may cur¬ 
rently overestimate the number of big meteoroids in the inter¬ 
planetary space, so they also may be improved by such a fit. 

Maps of the thermal emission of silicate particles are qualita¬ 
tively similar with those of the carbonaceous meteoroids and are 
not shown here. However, the emission spectra of the zodiacal 
cloud for the silicate particles are discussed in the next Section. 

We make maps of the thermal emission from the iron me¬ 
teoroids since they are remarkably different from the two ma¬ 
terials discussed above (Fig. M- Fine dust grains bulking near 
the Sun remain the brightest source regardless of the wavelength 
of observation. This stems from a very weak dependence of the 
absorption efficiency of iron particles on their size (within the 
range considered in this paper, see Fig. [9]). 

The difference in morphology of the maps for carbonaceous 
and iron particles suggests the way to discriminate between dif¬ 
ferent sorts of interplanetary meteoroids from the microwave ob¬ 
servations, when their accuracy permits. 
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Fig. 14. All-sky maps of the thermal emission from the interplanetary dust predicted by two meteoroid engineering models, assuming that their 
constituent particles are composed of iron. The reference frame is as in Fig.[T2] The brightness is provided in units of MJy sterad -1 for the COBE 
waveband centered at 240 pm, and in units of pK for the Planck and WMAP wavebands. 
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Fig. 15. Thermal emission spectra of the Zodiacal cloud predicted by 
the Kelsall model for the apex of Earth’s motion about the Sun and 
ecliptic pole. Large squares are the excess temperatures w.r.t. the black- 
body emission of the CMB (right scale), small squares are the absolute 
intensities of emission (left scale). 


4.2. Thermal Emission Spectra 

Figures [15] [16] and [IT] show the thermal emission spectra of 
the Zodiacal cloud predicted by Kelsall et al., pointed toward 
the apex of Earth’s motion about the Sun and ecliptic pole. 
The graphs are provided with two brightness scales: the abso¬ 
lute intensities are more natural and easier to read in the in¬ 
frared wavelengths, whereas the excess temperatures are instan¬ 
taneously comparable with other foreground emission sources 
in the CMB studies. The excess temperatures AT are defined 
in accord with iDikarev et al] (120091) so that B,i(Tcmb + A T) = 
Ba(Tcmb) + e, where e is the intensity of emission in excess of 
the CMB radiation at the temperature 7 cmb = 2.725 K, with 
the blackbody intensity B\ being calculated precisely with the 
Planck formula rather than in the Rayleigh-Jeans approximation. 


The Kelsall model is evaluated at the wavelengths of obser¬ 
vations with COBE DIRBE and Planck HFI, for which the ab¬ 
sorption efficiencies (“emissivity modification factors”) of the 
dust particles in the model components have been determined, 
and using the extrapolations described in Sect. l3.H for the WMAP 
wavelengths. The model predicts rather low emission in the mi¬ 
crowaves, e.g. below 1 pK in the apex and below 0.1 pK in the 
ecliptic poles for A > 3 pm. The emission in the apex is almost 
entirely due to the asteroid dust bands presumably having flat 
absorption efficiencies at these wavelengths, while a continuing 
decrease of the excess temperature with the wavelength increase 
for the polar line of sight for T > 3 mm indicates that the smooth 
background cloud is contributing with (Libs cx: A 2 . According 
to the Kelsall model, the foreground emission due to the Zodia¬ 
cal dust could indeed be disregarded in the WMAP survey of the 
microwave sky with its accuracy of 20 pK per 0.3° square pixel 
over one year. 

The Divine model in Fig. [16] turns out to be a lot brighter, 
however. It predicts an excess temperature of up to ~ 20 pK in 
WMAP’s W Band at A = 3 mm, and ~ 5 pK in the K Band 
at A — 13.6 mm, assuming that the meteoroids are all carbona¬ 
ceous. Silicate meteoroids would add much less than 10 pK in 
the entire WMAP wavelength range. If the interplanetary dust 
particles were composed exclusively of iron, the Divine model 
would not let them shine brighter than the Kelsall model. 

Predictions made with IMEM for the same compositions as 
above are shown in Fig. M IMEM is dimmer than the Divine 
model, it is still brighter than that by Kelsall. 

We have also estimated the total mass and cross-section of 
meteoroids in the Divine model and IMEM within 5.2 AU from 
the Sun, i.e. inside a sphere considered by iKelsall et al.l d 1998b 
as the outer boundary of the Zodiacal cloud. The Divine model 
contains ~ 2 • 10 2l) g mass and 6 • 10 6 AU 2 cross-section area 
of meteoroids. Assuming a density of 2.5 g cm -3 , the total mass 
of meteoroids could be contained in a single spherical body of 
the radius ~ 30 km, the total cross-section area correspons to a 
single sphere of the radius ~ 2 ■ 10 5 km. The last value of radius 
would be possessed by a planet three times bigger than Jupiter. 
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Fig. 16. Thermal emission spectra of the Zodiacal cloud predicted by 
the Divine model for the apex of Earth’s motion about the Sun and 
ecliptic pole, assuming that the dust particles are composed of amor¬ 
phous carbon (top), silicate (middle) and iron (bottom) material. Thick 
curves are the excess temperatures w.r.t. the blackbody emission of the 
CMB (right scale), dashed curves are the absolute intensities of emis¬ 
sion (left scale). 


IMEM has only 6 • 10 19 g mass and 6 • 10 6 AU 2 cross-section 
area of meteoroids within 5.2 AU from the Sun. 

iFixsen & Dwekl d2002h assumed a certain size distribution 
of interplanetary meteoroids and used the optical properties for 
silicate, graphite and amorphous carbonaceous particles in or¬ 
der to reproduce the annually-averaged spectrum of the Zodiacal 
cloud measured with the FIRAS instrument on board COBE. The 
smooth cloud of the Kelsall model was used to describe the spa¬ 
tial number density distribution of meteoroids. They found the 
total mass of meteoroids in the range 2-11 -10 1 s g, i.e. remark¬ 
ably lower than in the Divine model and IMEM. 


Fig. 17. Thermal emission spectra of the Zodiacal cloud predicted by 
IMEM for the apex of Earth’s motion about the Sun and ecliptic pole, 
assuming that the dust particles are composed of amorphous carbon 
(top), silicate (middle) and iron (bottom) material. Thick curves are the 
excess temperatures w.r.t. the blackbody emission of the CMB (right 
scale), dashed curves are the absolute intensities of emission (left scale). 


5. Conclusion 


We have made predictions of the microwave thermal emis¬ 
sion from the Zodiacal dust cloud in the wavelength range 
of COBE DIRBE, WMAP and Planck observations. We 
used the iKelsall et ali ( 199 80 model extrapolated to the mi¬ 
crowaves using the optical properties of interplanetary dust in¬ 
ferred bv lPlanck Collaboration XIVl (12014l) and lFixsen & D wek 
(200 2j), and two en gineering meteoroid models (IDivine 11993: 


Dikarev et al. 2004) in combination with the Mie light-scattering 
theory to simulate the thermal emission of silicate, carbonaceous 
and iron spherical dust particles. 

We have found that the meteoroid engineering models depict 
the thermal emission substantially brighter and distributed dif¬ 
ferently across the sky and wavelengths than the Kelsall model 
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does, due to the presence of large populations of macroscopic 
particles in the engineering meteroid models that are only avail¬ 
able in the asteroid dust bands of the Kelsall model. Both the Di¬ 
vine model and IMEM confirm an earlier es tima te of a ~ 10 pK. 
thermal emission from interplanetary dust dDikarev et al.ll2009k 
for the WMAP observations, provided that the dominant particle 
composition is carbonaceous. At smaller solar elongations, in¬ 
terplanetary dust can be orders of magnitude brighter, naturally. 

More detailed search for and account of interplanetary dust 
are therefore worthwhile in the CMB experiment results. Maps 
of the microwave sky should also be scrutinized to constrain the 
number of macroscopic particles in the engineering meteoroid 
models, as they are important to assess the impact hazard for 
long-term manned missions with large spaceships to the Moon 
or Mars. 
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